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Abstract 

A review of fundamentals and physical applications of fractional quan- 
tum mechanics has been presented. 

Fundamentals cover fractional Schrodinger equation, quantum Riesz 
fractional derivative, path integral approach to fractional quantum me- 
chanics, hermiticity of the Hamilton operator, parity conservation law and 
the current density. Applications of fractional quantum mechanics cover 
dynamics of a free particle, new representation for a free particle quantum 
mechanical kernel, infinite potential well, bound state in 5-potential well, 
linear potential, fractional Bohr atom and fractional oscillator. 

We also review fundamentals of the Levy path integral approach to 
fractional statistical mechanics. 

PACS number(s): 05.40.Fb, 05.30.-d, 03.65. Sq 



1 Introduction 

Classical mechanics and quantum mechanics are based on the assumption that 
the Hamilton function has the form 

tf(p,r) = fjUm (1) 

where p and r are the momentum and space coordinate of a particle with 
mass m, and V(r) is the potential energy. In quantum mechanics, p and r should 
be considered as quantum mechanical operators p and f . Then the Hamiltonian 
function H(p,r) becomes the Hamilton operator H(p,r), 

H( P: ?) = ^ + V(r), (2) 
where V(r) is the potential energy operator. 
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The square dependence on the momentum in Eqs.([T|) and ^ is empirical 
physical fact. However an attempt to get insight on the fundamentals behind 
this fact posts the question: are there other forms of kinematic term in Eqs.([T]) 
and @ which do not contradict the fundamental principles of classical me- 
chanics and quantum mechanics? A convenient theoretical physics approach to 
answer this question is the Feynman path integral approach to quantum me- 
chanics PP, as it was first observed by Laskin [2]. Indeed, the Feynman path 
integral is the integral over Brownian-like paths. Brownian motion is a special 
case of so-called a-stable probability distributions developed by Levy [3] and 
Khintchinc [4 . In mid 1930's they posed the question: Does the sum of N in- 
dependent identically distributed random quantities X = X\ -\-X-i... + Xn have 
the same probability distribution pn{X) (up to scale factor) as the individual 
steps pi(Xi), i = 1, ...AT? The traditional answer is that each pi(Xi) should be 
a Gaussian, because of the central limit theorem. In other words, a sum of N 
Gaussians is again a Gaussian, but with N times the variance of the original. 
Levy and Khintchine proved that there exist the possibility to generalize the 
central limit theorem. They discovered a class of non-Gaussian a-stable (stable 
under summation) probability distributions. Each a-stable distribution has a 
stability index a, often called the Levy index < a < 2. When a = 2 Levy 
motion is transformed into Brownian motion. 

An option to develop the path integral over Levy paths was discussed by Kac 
[5] , who pointed out that the Levy path integral generates a functional measure 
in the space of left (or right) continuos functions (paths) having only disconti- 
nuities of the first kind. The path integral over the Levy paths has first been 
introduced and elaborated with applications to fractional quantum mechanics 
and fractional statistical mechanics by Laskin (see, [2], [H]-[l2])- He followed 
the framework of the Feynman space-time vision of quantum mechanics, but 
instead of the Brownian-like quantum mechanical trajectories, Laskin used the 
Levy-like ones. If the fractal dimension (for definition of fractal dimension, see 
[S], 0) of the Brownian path is d^^^ m =2, then the Levy path has fractal 
dimension d|^^j = a, where a is so-called the Levy index, 1 < a < 2. The 
Levy index a becomes a new fundamental parameter in fractional quantum and 
classical mechanics similar to difT?!?™ 8 = 2 being a fundamental parameter in 
standard quantum and classical mechanics. The difference between the fractal 
dimensions of the Brownian and Levy paths leads to different physics. In fact, 
fractional quantum dynamics is generated by the Hamiltonian function H a (p, r) 
of the form g], [H]-[I7| 

H a (p,r) =D a \p\ a + V(r), Ka<2, (3) 

with substitutions p >p, r > r, and D a is the generalized coefficient, the 
physical dimension of which is [D a ] = erg 1 ^" • cm" • sec~ Q . One can say that 
Eq.j3|) is a natural generalization of the well-known Eq.([T]). When a = 2, 
D a = l/2m and Eq.([3]) is transformed into Eq.|T|) [2 . As a result, the fractional 
quantum mechanics based on the Levy path integral generalizes the standard 
quantum mechanics based on the well-known Feynman path integral. Indeed, if 
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the path integral over Brownian trajectories leads to the well-known Schrodinger 
equation, then the path integral over Levy trajectories leads to the fractional 
Schrodinger equation. The fractional Schrodinger equation is a new fundamental 
equation of quantum physics and it includes the space derivative of order a 
instead of the second (a = 2) order space derivative in the standard Schrodinger 
equation. Thus, the fractional Schrodinger equation is the fractional differential 
equation in accordance with modern terminology (see, for example, [9]-|13]). 
This is the main point for the term fractional Schrodinger equation or for more 
general term fractional quantum mechanics, FQM [2], |14j . When Levy index 
a = 2, Levy motion becomes Brownian motion. Thus, FQM includes standard 
QM as a particular Gaussian case at a = 2. The quantum mechanical path 
integral over the Levy paths [2] at a = 2 becomes the Feynman path integral 

In the limit case a = 2 the fundamental equations of fractional quantum 
mechanics are transformed into the well-known equations of standard quantum 
mechanics Q], [5], UM- 

2 Fractional Schrodinger equation 

2.1 Quantum Riesz fractional derivative 

Equation ([I]) lets us conclude that the energy E of a particle of mass m under 
the influence of the potential V(r) is given by 

E=^+V(r). (4) 

To obtain the Schrodinger equation we introduce the operators following the 
well-known procedure, 

d 

E ^ ih—, p ^ -ffiV, (5) 

where V = d/dr and h is the Planck's constant. Further, substituting 
transformation (JSJ into Eq.([T]) and applying them to the wave function i/j(r,t) 
yields 

ih ^dT- = "£ A ^ (r ' <} + v ^ r ^ ( 6 ) 

here A = V ■ V is the Laplacian. Thus, we obtain the Schrodinger equation 

i- 

By repeating the same consideration to Eq. (|3]) we find the fractional Schrodinger 
equation [2], [16] 

ih djA^^ D 2 /2 v Ka<2, (7) 

at 
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with 3D generalization of the fractional quantum Riesz derivative (— h, A) a > 2 
introduced by 

(-ft 3 A)«/V(r,t) = J dV ¥ |p|>(p,*), (8) 

where the wave functions in space ^>(r, £) and momentum <p(p,t) represen- 
tations are related each to other by the 3D Fourier transforms 



V>M) = ^3 J dVfrp(p,t), v(P,*) = _/ dVe-^frt). (9) 

The 3D fractional Schrodinger equation Eq.([7]) has the following operator 
form 

where fractional Hamilton operator H a (j>, r) results from Eq.Q with quantum- 
mechanical operators p and r substituted instead of p and r, 

H a (p,r) =D a \p\ a +V(r), 1< a < 2. (10) 
The ID fractional Schrodinger equation has the form [2], |14 ) -[l6 ] 

ih fy<P>t) = - D Ch^) a i/>(x, t ) + V(x)^(x, t), Ka<2, (11) 
at 

where (KS7) a is the quantum Riesz fractional derivativ^E The quantum 
Riesz fractional derivative is defined by the following way [2J, [T4] 

oo 

W(x,f)=-^ | dpe^lp^v&t), (12) 

— OO 

where ip(p,t) is the Fourier transform of the wave function ifj(x^t) given by 

oo 

tp(p,t)= J dxe~ l ^ip(x,t), (13) 

— OO 

and reciprocally 

oo 

ip(x,t) = — / dpe l ^(p(p,t). 

— oo 

It is easy to sec that Eq. (|11[) can be rewritten in the operator form, namely 



lr The Riesz fractional derivative was originally introduced in |23| 
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ih^ = H a ^j, (14) 
where H a is the fractional Hamilton operator 

H a = -D a (hV) a + V(x). (15) 

For the special case when a = 2 and Z?2 = l/2m (see, for details [2], [14]), 
where m is the particle mass, Eqs.Q and (fTT|) are transformed into the well- 
known 3D and ID Schrodinger equations [8]. 

2.2 The hermiticity of the fractional Hamilton operator 

The fractional Hamiltonian H a given by Eq. (|15[) is the Hermitian operator in 
the space with scalar product 

00 

(<t>,X)= J dxct)*(x,t)x(x,t). (16) 
—00 

To prove the hermiticity of H a let us note that in accordance with the 
definition of the quantum Riesz fractional derivative given by Eq. (|12p there 
exists the integration-by parts formula 

(<t>,(hV) a x) = ((hV) a <f>, x ). (17) 

The average energy of a fractional quantum system with Hamiltonian H a is 



E a = J dxip*(x,t)H a <ip(x,t). (18) 

— OO 

Taking into account Eq. ffTTj) we have 

E a = I dx-4>*{x,t)H a ^[x,t)= I dx(H+tlj(x,t))*ip(x,t) = E*. 



As a physical consequence, the energy of the system is real. Thus, the 
fractional Hamiltonian H a defined by Eq. lfTS"]) is the Hermitian or self-adjoint 
operator in the space with the scalar product defined by Eq. (jT6)) [15], [16] 



(£#&X) = (0,tfaX). (19) 

The generalization of the proof of hermiticity for 3D case is straightforward. 
Note that Eq. fTTj) leads to the important equation 

^ I dxrP* (x,W(x, t) =0, (20) 
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which shows that the wave function remains normalized, if it is normalized 
once. Multiplying Eq. flTTt from the left by tp*(x,t) and the conjugate complex 
of Eq. fTTj) by ip(x,t) and then subtracting the two resultant equations finally 
yield 

d 

ih— (tp (x,t)il)(x,t)) = tp {x,t)H a <ip{x,t) - i/j(x,t)H a i/j (x,t). 

Integrating this relation over all values of the space variable and using the 
fact that the operator H a is self-adjoint, we find Eq. ([20|) . The above consider- 
ation can be easily generalized to 3D case. 

2.3 The parity conservation law for the fractional quan- 
tum mechanics 

It follows from the definition (|12p of the quantum Riesz fractional derivative 
that 

(ftV) Q exp{z^} - -\p\<* exp^}, (21) 

in other words, the function exp{ipx/h} is the eigenf unction of the quantum 
Riesz fractional operator (7iV) Q with eigenvalue — \p\ a . 
The 3D generalization is straightforward, 

(-h 2 A)^ 2 exp{z^} = |p| Q exp{z^}, (22) 

the function exp{ipx/?i} is the eigenfunction of 3D quantum Riesz fractional 
operator (— h 2 A) a / 2 with eigenvalue |p| Q . 

Thus, the operators (?iV) Q and (— h 2 A) a / 2 are the symmetrized fractional 
derivative, that is 

(hV x ) a ... = (W_ x ) a ..., (23) 

(-h 2 A r ) a / 2 ... = (-h 2 A_ r ) a / 2 .... (24) 

Because of the properties (|17p and P^)l the fractional Hamiltonian H a (see, 
for example Eqs.([5]) or ([5])) remains invariant under inversion transformation. 
Inversion, or to be precise, spatial inversion consists in the simultaneous change 
in sign of all three spatial coordinates 

r — > r, x — > —x, y —> —y, z — > —z. (25) 

Let us denote the inversion operator by P. The inverse symmetry is the fact 
that P and the fractional Hamiltonian H a commute, 

PH a = H a P. (26) 

We can divide the wave functions of quantum mechanical states with a well- 
defined eigenvalue of the operator P into two classes; (i) functions which are 



G 



not changed when acted upon by the inversion operator, Pip + (r) — "0+( r ) the 
corresponding states are called even states; (ii) functions which change sign 
under the action of the inversion operator, P , 0_(r) = —ip_(r) the corresponding 
states are called odd states. Eq. ([2"6l express the "parity conservation law" for 
the FQM [T6] ; if the state of a closed fractional quantum mechanical system has 
a given parity (i.e. if it is even, or odd) , then this parity is conserved. 

2.4 The current density 

By multiplying Eq.([7]) from left by ip*(r, t) and the conjugate complex of Eq.([7]) 
by ip{v,t) and subtracting the two resultant equations yield 

/>r^*(r,t)^(r,*)) = (27) 



^ f d\ (^(r,i)(-?i 2 Ar/V(r,^-VXr,t)(-n 2 A)«/V(r,t)) 



From this integral relationship we are led to the following well-known differ- 
ential equation 

^M+diyj(r,i) = 0, (28) 

where p(r, t) — ip* (r, t)ip(r, t) is the quantum mechanical probability density 
and the vector j(r, t) can be called by the fractional probability current density 
vector 



j(r,t) = — (7/.*(r,i)(-?i 2 A) a / 2 - 1 V^(r,t) - ^(r, t)(-?i 2 A)"/ 2 " 1 Vi/>*(r, i)) , 

(29) 

where we use the following notation V = d/dr. Introducing the momentum 
operator p = j ^ we can write the vector j in the form [16] 

j =D a (^(p 2 )"/ 2 - 1 ^* + r(p* 2 ) a/2 - 1: P^) (30) 

The new fundamental Eqs. (|2"9"]) and (T5UI) are the fractional generalization 
of the well-known equations for probability current density vector of standard 
quantum mechanics [5]. 

To this end, we express Eq. (|5(I)) in the terms of the velocity operator, which 
is defined as follows v = dr/dt, where r is the operator of coordinate. Using the 
general quantum mechanical rule for differentiation of operator 




we have 
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v = -r( H ar - rH a ), 
n 

Further, with help of the equation f(p)r — rf(p) = —ihdi/dp, which holds for 
any function f(p) of the momentum operator, and taking into account Eg. (1101) 
for the Hamiltonian operator H a (p,r) we obtain the equation for the velocity 
operator 

v = aD a \p*\*/ 2 - 1 p, (31) 

here p is the momentum operator. By comparing of Eqs. (|3"0")) and (|3"Tj) we 
finally conclude that 

j =- (ipvtp* + ip*v^p) , 1 < a < 2. (32) 
a 

To get the probability current density equal to 1 (the current when one 
particle passes through unit area per unit time) the wave function of a free 
particle has to be normalized as 

<P(r,t) = ^exp{^pr-^Et}, E = D a \p\ a , Ka<2, (33) 

where v is the particle velocity, v = aD a p Q_1 . Then we have 

j=~, v = aD a \p 2 \%- 1 p, (34) 
that is, the vector j is indeed the unit vector. 

Equations (|2T))) - (f3"4")) are the fractional generalization of the well-known equa- 
tions for probability current density vector and velocity vector of the standard 
quantum mechanics [8]. 

2.5 The time-independent fractional Schrodinger equation 

The special case when the Hamiltonian H a does not depend explicitly on the 
time is of great importance for physical applications. It is easy to see that in 
this case there exist the special solution of the fractional Schrodinger equation 
HU of the form 

i>{x,t) =e-W Et cj>{x), (35) 

where <j)(x) satisfies 

H a (f>(x) = Ecj){x), (36) 

or 

~ D a {hV) a <j){x) + V{x)<j){x) =Ecf>(x), l<a<2. (37) 
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The equation ([3"T]) we call by the time-independent (or stationary) fractional 
Schrodinger equation [15], [16]. We see from Eg . ([35]) that the wave function 
4>{x,t) oscillates with a definite frequency. The frequency with which a wave 
function oscillates corresponds to the energy. Therefore, we say that when the 
fractional wave function i})(x,t) is of this special form, the state has a definite 
energy E. The probability to find a particle at x is the absolute square of the 
wave function \ip\ 2 . In view of Eq. (l35]) this is equal to |0| 2 and does not depend 
upon the time. That is, the probability of finding the particle in any location is 
independent of the time. We say under these circumstances that the system is 
in a stationary state - stationary in the sense that there is no variation in the 
probabilities as a function of time. 



3 Path integral 

If a particle at an initial time t a starts from the point x a and goes to a final 
point Xb at time tb, we will say simply that the particle goes from a to b and its 
trajectory (pathj^l x(t) will have the property that x(t a ) = x a and x(tb) = Xb- 
In quantum mechanics, then, we will have an quantum-mechanical amplitude, 
often called a kernel, which we may write K(xbtb\x a t a ), to get from the point a 
to the point b. This will be the sum over all of the trajectories that go between 
that end points and of a contribution from each pQ. For the one dimensional 
version of the Hamiltonian Eq. © 

H a {p,x) = D a \p\ a + V{x,t), (38) 

following consideration provided in [2], |14j . we come to the definition of the 
kernel K(xbtb\x a t a ) in terms of path integral in the phase space representation 

oo oo 

K(x b tb\x a t a ) = lim / dxi...dx N -i — / dpi...dp N x (39) 

N-too J (ZTtn) J 



\ i N \ i N i N 

exp | frX^'fo ~ x i- 1 ) \ x ex P j ~ft Da£ Yl - ^X^^J'e) 

here e = (tb — t a )/N , Xj = x(js), pj = p(js) and Xq = x a , xn = Xb- Then 
in the continuum limit N — > oo, e->0we have 

K L {x b t b \x a t a ) = (40) 



x(t b )=x b 



J Bx(t) J Dp(r)exp j % - J dr[p(T)x(r) - H a (p(r), x(r), r)] 



x(t a )=X a \ t a 

2 For simplicity, here we consider one dimensional motion. 
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where x denotes the time derivative, H a (jp(r), x(t), t) is the fractional Hamil- 
tonian given by Eq . with the replacement p — >p(r), x — > x(r) and {p(t),x(t)} 
is the particle trajectory in phase space, and, finally, the phase space path intc- 

x(ti,)=x b 

gral J Dx(t) J Dp(r)... is given by 

x{t a )—x a 

x(t b )=x b 



J Dx{t) J Dp(r)... = (41) 



c(t a )—x a 



lira / dxt...dx N -\ \ N I dpi...dpN* 



(.pi(xi-x a ) .D a e\pi\ a \ ( .p N (x b - x N -i) .D a e\p N \ a \ 
exp V h »— jj— jx...xexp|, 

The exponential in Eq. (|4"0"]) can be written as exp{j-S a (p, x)} if we introduce 
the fractional canonical classical mechanical action S a (p,x) for the trajectory 
p(t), x(t) in phase space [14] 

h 

S a (p,x) = J dr(p(r)x(T)-H a (p(T),x(T),T). (42) 

ta 

Since the coordinates xq, xn in the definition (|4ip are fixed at their initial 
and final points, £o = x a and xn = Xb, the all possible trajectories in Eq. (|40|) 
satisfy the boundary condition x(tb) — x b , x(t a ) = x a . We see that the definition 
given by Eq. (|41l) includes one more p^-integrals than Xj-integrals. Indeed, while 
xq and xjv are held fixed and the ^'-integrals are done for j = 1,...,N — 1, 
each increment Xj — Xj-% is accompanied by one -integral for j = 1, ...,JV. 
The above observed asymmetry is a consequence of the particular boundary 
condition. Namely, the end points are fixed in the position (coordinate) space. 
There exist the possibility of proceeding in a conjugate way keeping the initial 
p a and final p momenta and fixed. The associated kernel can be derived going 
through the same steps as before but working in the momentum representation 
(see, for example, [24]). 

The kernel K{xbtb\x a t a ) introduced by Eq. ([40]) describes the evolution of 
the quantum mechanical system 

oo 

ijj{xb,t b ) = J dx a K(xbtb\x a t a )ip(x a ,t a ), (43) 

— oo 

where ip(x a ,t a ) is the wave function of the initial state (at t = t a ) and 
tp(xb, tb) is the wave function of the final state (at t = tb). 
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3.1 Free particle 



For a free particle when V(x,t), we have H a (p) = D a \p\ a and it's easy to see 
that Eq.® results in 



K^(x b t b \x a t a ) = ^ [dp.e^\i P -^^-i D ^ b - ta) \. (44) 



here (x b t b \x a t a ) stands for a free particle kernel. 

Taking into account Eq.(j3U) it is easily to check on directly the consistency 
condition 



K^\x b t b \x a t a ) = / dx'K^\x b t b \x't')-K^(x't'\x a t a ) 



This is a special case of the general quantum-mechanical rule: for events 
occurring in succession in time the amplitudes are multiplied 

oo 

K{x b t b \x a t a ) = J dx'K{x b t b \x't')-K(x't'\x a t a ). (45) 



3.1.1 Fox ii-function representation for a free particle kernel 

Let's show how a free particle fractional quantum mechanical kernel (xbtb\x a t a ) 
defined by Eq. ([44"l) can be expressed in the terms of the Fox H- function [18], 
[19], [20]. Follow by [17], we obtain the Mellin transform of the quantum me- 
chanical fractional kernel defined by Eq. (|44)) . Comparison of the inverse Mellin 
transform with the definition of the Fox function yields the desired expression 
in terms of "known" function, i.e. Fox H- function. Note that H- function bears 
the name of its discoverer Fox [18] although it have been known at least since 
1888, according to [T9] . 

Introducing for simplicity the notations x = x b — x a , r = t b — t a , we rewrite 
Eq.flU]) 

OO 

^<* T > - 55 /*■•"{'£ "'^J- (46) 

— oo 

It is easy to see that the relation K^(x,t) = K^°'{—x,t) holds. Hence, 
it is sufficient to consider K^\x,t) for x > only. Further, we will use the 
following definitions of the Mellin transform 

A °° 

K( 0) (s,t) = J dxx s - 1 K { - Q \x,T), (47) 
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and inverse Mellin transform 

C+20C 



K^(x,t) = —. [ dsx~ s K^{s,T), (48) 
2m J 



c—ioo 



where the integration path is the straight line from c — zoo to c + ioo with 
< c < 1. 

The Mellin transform of the K^°'(x,t) defined in accordance with Eq. (|4l)]) 

is 

oo oo 



By changing of the variables of integration p — >• I jj^-p j c, x — > ( jD r J 
one obtains the integrals in the complex <r and £ planes. Considering the paths 

A 

of integration in the c and £ planes it is easy to represent K^°\s,t) as follows 

A 



OO OO 

J df£ s - x J dcexp{^-|c| a }- (49) 

-oo 



2tt \{h/iD a T) x / a , 
The integrals over <i£ and <i<r can be evaluated by using the equation 

oo oo 

J die^ J dc exp{^ -c a } = 7 4y sin ^^r( S )r(l-^), (50) 



where s — 1 < a < 2 and T(s) is the gamma function^. 
Inserting Eq. (l50|) into Eq. (|49|) and using the functional relations for the 
gamma function, T(l — z) = —zT(—z) and T(z)r(l — z) — tt/ sin7T2;, yield 

J&(.,r) 1 ' » V "' r<S)r( 



The inverse Mellin transform gives a free particle quantum mechanical kernel 



C+200 



= | dsx- s K^( S ,r) 



3 The gamma function T(s) has familiar integral representation T(s) = J dtt s 1 e *, Res > 



0. 
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2nia J \{h/iD a Ty/ a ) r(±^)r(±±s)' 

c— zoo 

where the integration path is the straight line from c — ioo to c + ioo with 
< c < 1. Replacing s by — s we obtain 



(51) 



a \ {'h/iD a T) 1 / a J 2m 



-c+ioc 



ds 



1 



h \iD a T 



1/a 



r(-*)r(^) 
r(i±a)r(i=2) 



— c— zoo 



The path of integration may be deformed into one running clockwise around 
R + — c. Comparison with the definition of the Fox H- function (see, Eqs.(58) 
and (59), in [T7])) leads to 



A-<°>(aV 



(52) 



a \(H/iD a r) 1 / c 



H. 



i.i 



1 



ft, \iD a T 



1/a 



(1- 1/a, 1/a), (1/2, 1/2) 
(0,1), (1/2, 1/2) 



Applying the Property 12.2.5, Ref.|17]. of the Fox ff-function we can express 
(x,t) as 



1 ( h V /a (1, 1/a), (1,1/2) 
ft\iD aT ) 1 (1,1), (1,1/2) 



x > 0. (53) 



i.i 



Or for any x, 

kW(x,t) = —II 
a 

Substituting x = x b — x ai t = t b — t a , finally yields 

K^(x b t b \x a t a ) = 



If ft \ 1/a (1, 1/a), (1,1/2) 
' " (LI), (1,1/2) 



ft \lD a T 



(54) 



a\Xb - x a \ 



i.i 



iD a (t b - t a ) 



1/a 



\Xb - X a \ | 



(1, 1/a), (1,1/2) 
(1,1), (1,1/2) 



(55) 



This is new equation for ID free particle fractional quantum mechanical 
kernel K ^ (x b t b \x a t a ). 
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Let us show that Eq. ([55j) includes as a particular case at a = 2 the well- 
known Feynman quantum mechanical kernel, see Eq.(3-3) in [1 ]. Setting in 
Eq. (|55l) a = 2, applying the series expansion for the function 

^ 2 [n \iD a (t b -t a )J ]Xb Xal] (i,i),(i,i/2) _ ' 

and finally, substituting k — > 21 yield 



1/2 \ 2i W 21 



(56) 



Taking into account the identity T(^ + z)T(^ — z) — co J 7rz , and applying the 



Gauss multiplication formula r(2i) = y^o — r(Z)r(Z + i), we find that 



2tt 



(20!r(- -0 = ^(2)^!. (57) 
With help of Eq. (f57)) the kernel K^°\x, r)\ a =2 can be rewritten as 



r)i 2 - - j- m 1/2 f f-i f±yT H)l|f 



I / ft. \ 1/2 1 Id 2 

1 exp { -i#-}. 



(58) 



2y/7rh \iD 2 r J AfiiD 2 T' 
Since £>2 = l/2m we come to the Feynman kernel (see Eq.(3-3), PQ) 

Thus, it is shown how Feynman a free particle kernel can be derived from 
the general equation (|55j) for the fractional quantum mechanical kernel. 



4 Applications of fractional quantum mechanics 

4.1 A free particle fractional Schrodinger equation 

4.1.1 Scaling properties of ID fractional Schrodinger equation for a 
free particle 

To make general conclusions regarding solutions of ID fractional Schrodinger 
equation for a free particle, let's study the scaling of the solutions. The scale 
transformations could be written as 
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t = At', x = AV, D a = y D' a , ijj{x, t) = \ 5 ^{x', t'), 



where /3, 7, S are exponents of the scale transformations which should leave 
invariant a free particle ID fractional Schrodinger equation 

in di^t)_ = _ Da{m)a%j}{x ^ (59) 

00 

and save the normalization condition J dx\ip(x,t)\ 2 — 1. It reduces the 

—00 

number of exponents up to 2 and brings the two-parameters scale transformation 
group 



t = Xt', x = \ f3 x', D a = \ af} - 1 D' al ij}{\ p x,H) = X~ 0/2 ijj(x,t), 

(60) 

where /3 and A are arbitrary group parameters. 

When the initial condition ijj(x,t — 0) is invariant under the scaling group 
Eq. (|60|) then the solution of Eq. (|59j) remains the group invariant. As an example 
of invariant initial condition one may keep in mind ip(x,t = 0) = S(x), which 
gives us the Green function of ID fractional Schrodinger equation. 

To get the general scale invariant solutions of ID fractional Schrodinger 
equation we may use the renormalization group framework. As far as the scale 
invariant solutions of Eq. ([59| should satisfy the identity Eq. (|60l) for any arbi- 
trary parameters /? and A, the solutions can depend on combination of x and 
t to provide the independency of /3 and A. Thus, because of existence of the 
relationships between scaling exponents, a(3 — 7 — 1 =0 and <5 + (3/2 = the 
solutions are 

^(x,t) = U(x/(D a t)i) = —}—^( x /( Da t)i), (61) 

X {-L>at) a 

where arbitrary functions $ and ^ are determined by the conditions, <!>(.) = 
V(l,.) andtf(.)=V(-,l)- 

4.1.2 Exact solution 

Following [15], [22] let's solve ID fractional Schrodinger equation for a free 
particle (|59j) with some initial condition ip^ix) 

iP(x,t = Q)=xP Q (x). (62) 

Applying the Fourier transforms Eqs. pU]) and using the quantum Riesz frac- 
tional derivative Eq. (|12[) yield for the wave function ip(p,t) in the momentum 
representation, 

m^M=B a b|^(p,t), (63) 
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with the initial condition (fo(p) given by 

oo 

<p o (p) = t p(p,t = 0)= [ dxe-^%(x). (64) 



The solution of the problem Eqs. (R)5]) and ([Ml) is 

tp(p,t) = exp{-i Da ^ t }ip (p), (65) 

Therefore, the solution of ID fractional Schrodinger equation Eq. ([S"9"]) with 
initial condition given by Eq. (|62[) can be presented as 



= 2^ / **> J dpexp{i P{x % X>) - iP^ly^l (66) 

— oo — oo 



r/>(x,t)= (67) 



oo 

/d,X -, 7\H 2 2 
a(x — x') 

— oo 

Here we expressed the integral over dp in Eq. (|B5| in terms of ilj '2 -function. 
If we choose the initial condition V'o( a; ) = ^o(^), then Eq. ([67]) gives us quan- 
tum mechanical kernel K^°\x,t\0,0) for ID free particle fractional Schrodinger 
equation 



iD a t 



l/a 



(1,1/a), (1,1/2) 
(1,1), (1,1/2) 



%(x'). 



K (0 \x,t\0,0) 



-H. 



i.i 



1 



h \iD a t 



l/a 



(1,1/a), (1,1/2) 
(1,1), (1,1/2) 



(68) 



or applying the Property 12.2.5 of the Fox iJ-function (see, [17]) we can 
write for (x,t\0,0) 



K^(x,t\0,0) 



(69) 



1 

a 



Xh/iD a r)^ 



H. 



i.i 



1 



h \iD a T 



l/a 



(1- l/a, l/a), (1/2, 1/2) 
(0,1), (1/2, 1/2) 



It is easy to see that Eqs. ([66l) and ([69]) are scale invariant solutions of ID 
fractional Schrodinger equation (see, Eq. ([59]) ) for a free particle. 
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4.1.3 3D generalization 

A free particle quantum dynamics in 3D is governed by following equation (see, 
Eq.©) 

iU ^Bt^ = D a (-H 2 A) a ^(r, t), ^(r, t = 0) = ^„(r). (70) 

Using 3D Fourier transformed defined by Eq.© and the definition of 3D 
quantum fractional Riesz derivative given by Eq.® yield for the wave function 
ip(p,t) in the momentum representation, 

in^^- = D a \p\^( P ,t), (7i) 

with the initial condition (p (p) given by 

oo 

yo(p) = V(P.* = 0)= / d 3 re-^^ (r). (72) 

— OO 

Go back to Eq. (fT0[) we can see that the solution ^(r, t) has a form 



— oo — oo 

The integral over d 3 p can be expressed in terms of H 3 '3 -function, see, for 
instance Eqs.(33) and (34) in.[T7]- Thus, the solution of the problem Eg. (1701 is 

#•,*)= (73) 



1 

2ira 



d 3 r'-, 



H. 



1.2 



r /|3 3.3 



1 



fi \ iD a t 



l/a 



r — r 



„ , (1,1), (1, l/a), (1,1/2) 
(1,1), (1,1/2), (2,1) 



^o(r')- 



Substituting into Eq. (|T3"|) ^ ( r ) = ^o( r ), gives us quantum mechanical kernel 
K^(r,t\Q,0) for a free particle 3D fractional Schrodinger equation 



Jf (o) (r,i|0,0) 



1 1 
2na |r 



TJ L ' Z 
1 H 3,3 



iD a t 



l/a 



(1,1), (1, l/a), (1,1/2) 
(1,1), (1,1/2), (2,1) 



(74) 

This is new equation for a free particle quantum mechanical 3D kernel. We 
see that in comparison with ID case 3D quantum kernel is expressed in the 
terms of H 3 ' 3 Fox H- function. In the case a — 2 we come to the well-known 

equation for Feynman 3D quantum kernel Kp\r, t\0, 0), 
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^^^.ifhM.^f^fy:}. ,75, 

4.2 The infinite potential well 

A particle in a one-dimensional well moves in a potential field V(x) which is 
zero for —a<x<a and which is infinite elsewhere, 

V(x) = oo, x < —a (i) 

V(x) =0, -a < x < a (ii) (76) 

V(x) = oo, a; > a (iii) 

It is evident a priori that the spectrum will be discrete. We are inter- 
ested in the solutions of the fractional Schrodinger equation ([37]) that describe 
the stationary states with well-defined energies. Such a stationary state with 
an energy E is described by a wave function if>(x,t) which can be written as 
tjj(x,t) = exp{— z-^}0(cc),where <p(x) is now time independent. In the regions 
(i) and (iii), (see Eq. ([76j) ) we have to substitute oo for V(x) into Eq. (j37|) . and it 
is easily to see that here the fractional Schrodinger equation can be satisfied only 
if we take 4>{x) = 0. In the middle region, (ii), the time-independent fractional 
Schrodinger equation is 

- D a {TiV) a <t>{x) = E<t>{x). (77) 

We can treat this as a fractional eigenvalue problem [15]. Within region (ii), 
the eigenfunctions are determined by Eq. (|77|) . Outside of the region (ii), x < — a 
and x > a, the eigenfunctions are zero. We want the wave function <fi(x) to be 
continuous everywhere, and this means that we impose the boundary conditions 
(/>(— a) = 4>{a) = for the solutions of the fractional differential equation (177)) . 
Then the solution of Eq.(|77|) in the region (ii) can be written as 

cj) cvcn (x) = Acoskx, or odd (x) = Asinkx, 
where the following notation is introduced 

k =hih x ' a > 1<a£2 - {78) 

The even (under the reflection x — > —x) solution 4> cvcn (x) satisfies the bound- 
ary conditions if 

fc=(2m+l)^, m = 0,1, 2,3,... 
2a 
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The odd (under the reflection x — » — x) solution ^> odd (a;) satisfies the bound- 
ary conditions if 



ft I -7T 

k= , m =1,2,3,... (79) 

a 

It is easy to check that the normalized solutions are 

9C en O) = -7= cos ! ( m + ib— 

ya L 2 a 

and 

/odd/ \ 1 • TO7r £ 

9 m i x ) = -7= sin . 

The solutions (/) cvcn (a;) and </> odd (x) have the property that 

a 



dxC en (^)C dd W = 



where S mn is the Kronecker symbol and 

The eigenvalues of the particle in a well with help of Eas.ff8| and (|79)l are 



ttTi 



n a , n= 1,2,3...., 1 < a < 2. (80) 



It is obviously that in the Gaussian case (a = 2) Eq. ([80|) is transformed to 
the standard quantum mechanical equation (for example, see Eq.(20.7), Ref.[8]) 
for the energy levels for a particle in a box. 

The state of the lowest energy, the ground state, in the infinite potential well 
is represented by the (j)^ en (x) at m — 0, 



1 r ttx , 



Aground fa) = C™^) = — C0S {^} 



and its energy is 



Aground = D a ( — J . (81) 
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4.3 Fractional Bohr atom 



When V(r) is the hydrogenlike potential energy 

where e is the electron charge, Ze is the nuclear charge of the hydrogenlike 
atom, we come to the eigenvalue problem for fractional hydrogenlike atom, 

D a (-h 2 Ar/^(r) - -^0(r) = E<f>(r). 

This eigenvalue problem had been solved at first in [15] ■ The total energy 
is E = Ekin + V, where Ekin is the kinetic energy Ekin — -D Q |p| Q , and V is 
the potential energy V = It is well-known that if the potential energy 

is a homogeneous function of the coordinates and the motion takes place in a 
finite region of space, there exists a simple relation between the time average 
values of the kinetic and potential energies, known as the virial theorem (see, 
page 23, [27]). It follows from the virial theorem that between average kinetic 
energy and average potential energy of the system with Hamiltonian ([3]) there 
exist the relation 



aE kin = -V, (82) 
where the average value / of any function of time is defined as 

oo 

/= lim i / dtf(t). 

T-s-oo 1 J 


In order to evaluate the energy spectrum of the fractional hydrogenlike atom 
let us remind the Niels Bohr postulates [25] , [2H] : 

1 . The electron moves in orbits restricted by the requirement that the angu- 
lar momentum be an integral multiple of h, that is, for circular orbits of radius 
a n , the electron momentum is restricted by 



pa n = nh, (n= 1,2,3, ...), (83) 

and furthermore the electrons in these orbits do not radiate in spite of their 
acceleration. They were said to be in stationary states. 

2. Electrons can make discontinuous transitions from one allowed orbit cor- 
responding to n — ni to another corresponding to n = ri2, and the change in 
energy will appear as radiation with frequency 



E — E 



(84) 



h 

An atom may absorb radiation by having its electrons make a transition to 
a higher energy orbit. 
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Using the first Bohr's postulate and Eq.(|8"2"j) yields aD a fa^ = from 

which it follows the equation for the radius of the fractional Bohr orbits 

On = a n~ , (85) 

here ao is the fractional Bohr radius (the radius of the lowest, n — 1, Bohr 
orbit) defined as, 

f aD a h a \^ 

fl o = ■ ( 86 ) 



V Ze 

By using Eq. (|52"|) we find for the total average energy, E = (1 — a)Eki n - 
Thus, for the energy levels of the fractional hydrogen-like atom we have 

E n = (1 - a)E Q n~^, 1< a < 2, (87) 

where Eq is the binding energy of the electron in the lowest Bohr orbit, that 
is, the energy required to put it in a state with E = corresponding to n = oo, 



a a DJi a , 

The energy (a — l)Eo can be considered as a generalization of the Rydberg 
constant of standard quantum mechanics. It is easy to see that at a = 2 the 
energy (a — 1)Eq is transformed into the well-known expression for the Rydberg 
constant, Ry = me i /2h 2 . 

The frequency of the radiation oj associated with the transition, say, for 
example from m to n, m — > n, is, according to the second Bohr postulate, 



(1 - a)E 
OJ — 



1 1 



(89) 



The new equations (I551 - (|89l) bring us fractional generalization of the "Bohr 
atom" theory. In the special Gaussian case (standard quantum mechanics) 
Eqs. (l55|) - (|M| allow us to reproduce the well-known results of the Bohr theory 



4.4 Fractional oscillator 

4.4.1 Quarkonium and fractional oscillator 

As an another physical application of the developed FQM we propose a new frac- 
tional approach to study the quark-antiquark qq bound states treated within the 
non-relativistic potential picture [2]. Note, that only for heavy quark systems 
(for example, charmonium cc or bottonium bb) the non-relativistic approach 
can be justified. The term quarkonium is used to denote any qq bound state 
system [28 in analogy to positronium in the e + c~ system. The non-relativistic 
potential approach remains the most successful and simplest way to calculate 
and predict energy levels and decay rates. 
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Thus, from stand point of " potential" view, we can assume that the confining 
potential energy of two quarks localized say, at the space points and Yj is given 

by 

V{\T t -T j \)=q i q j \T i -T j \P, (90) 

where qi and qj are the color charges of i and j quarks respectively and 
the index j3 > 0. Equation. (j90|) coincides with the QCD requirements: (i) At 
short distances the quarks and gluons appear to be weakly coupled; (ii) At large 
distances the effective coupling becomes strong, resulting in the phenomena of 
quark confinement^. 

Considering the TV-quarks statistical system yields the following equation for 
the potential energy U(r\, rjv) of the system, 

U{r 1 ,...,r N )= Yl (91) 

l<i<j<N 

In order to illustrate the main idea, we consider the simplest case, when color 

JV 

qi charge only can be q or — q and the colorless condition ^2 qi = takes place. 

i=i 

Using the general statistical mechanics approach (see the Definition 3.2.1 and 
the Proposition 3.2.2, Ref. [3T]) we can conclude that only for < (3 < 2 the 
many particle system with the potential energy (|9ip will be thermodynamically 
stable. 

In order to study the problem of quarkonium it seems reasonable to consider 
the non-relativistic FQM model with the fractional Hamiltonian operator H a ^ 
defined as 

Ha,p = D a (-Ti 2 A) a / 2 + g 2 |r|^, 1< a < 2, 1< < 2, (92) 

where r is 3D vector, A = d 2 /dr 2 is the Laplacian, and the operator 
{-h 2 A) a / 2 is defined by Eq.©. 

For the special case, when a = j3 the Hamiltonian operator (|92p has a form 

H a = D a {-h 2 A) a / 2 + q 2 \r\ a , l<a<2. (93) 

It is easy to see that the Hamiltonian H a is the fractional generalization of 
3D harmonic oscillator Hamiltonian of standard quantum mechanics. Follow [2], 
[16] . we will call by the fractional oscillators the quantum mechanical models 
with the Hamiltonians given by the Eqs. (|9"2l and (f9"3]) . 

It would be interesting to calculate the energies of the bound states and the 
decay rates based on the FQM models with the Hamiltonians and (j9"3"|) . 
These results could be compared with the experimental statistical data of J/ ip 
decays. The information on decay rates and angular distribution provides an 
ideal testing ground for the fractional models of qq bound states. 

4 The term quark confinement describes the observation that quarks do not occur isolated 
in nature, but only in hadronic bound states as the colorless objects such as baryons and 
mesons. 
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4.4.2 Spectrum of ID fractional oscillator in semiclassical approxi- 
mation 

The ID fractional oscillator with the Hamilton function H a — D a \p\ a + g 2 |a;|' 3 
poses an interesting problem for semiclassical treatment [TB]. We set the total 
energy equal to E, so that 

E = D a \p\ a + q 2 \xf, (94) 

(\ 1/q 
■j^-(E — q 2 \x\P)\ . At the turning points p = 0. Thus, 

classical motion is thus possible in the range \x\ < (E/q 2 ) 1 ^ 13 . 

A routine use of the Bohr-Sommerfeld quantization rule [5] yields 



2irh(n + -)= <j)pdx = 4 I pdx = 4 J D-^ a (E - q 2 \xf) la dx, (95) 

o o 

where the notation <f means the integral over one complete period of the 
classical motion, x m = (E/q 2 ) 1 ^ is the turning point of classical motion. To 
evaluate the integral in the right hand of Eq. (|95p we introduce a new variable 
y = x(E/q 2 )^ 1 / 13 . Then we have 



J D-^(E-q 2 \xfrdx = - T ±^EiH jdy(l-y^. 
o a q 

The integral over dy can be expressed in the terms of the B-function. Indeed, 
substitution z = y@ yields^] 

l l 
J dy(l - /)V« = I J dzz^tl - z)i = ~ + 1). (97) 

o o 
With help of Eq.^FQ we rewrite Eq.([|55]) as 

2T:h(n+-) = — E^+^-B{-,- + \). 

The above equation gives the values of the energy of stationary states for 
ID fractional oscillator [TBI. 



5 The B-function is defined by 

1 

i{u,v) = f dyy^il-yy- 1 . (96) 



B( 
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This equation is generalized the well- known energy spectrum of the standard 
quantum mechanical oscillator (see for example, [8]) and is transformed to it at 
a = 2, 13 = 2. 

We note that at 



Eq. (j551) gives the equidistant energy spectrum. When 1 < a < 2 and 1 < (3 < 2 
the condition given by Eq. (|99l) takes place for a — 2 and (3 = 2 only. It 
means that only standard quantum mechanical oscillator has equidistant energy 
spectrum. 

5 Some solvable models of fractional quantum 
mechanics 

Let's outlook a few analytically solvable problems of fractional quantum me- 
chanics recently elaborated in [25] . 

5.1 Bound state in 5— potential well 

For one dimensional attractive ^-potential well, V(x) = —jS(x), (7 > 0), where 
S(x) is the Dirac delta function, fractional Schrodinger equation Eq. ([37]l becomes 



In any one-dimensional attractive potential there will be a bound state, that 
is E < 0. In this case Dong and Xu [25J, found the energy and the wave function 
of the bound state. The bound energy has a form 



D a {hV) a <f)(x) - yS(x)(t>(x) = E<t>{x), 



1 < a < 2. 



(100) 




1 < a < 2, 



(101) 



here B(l/a, 1 — 1/a) is the i3-function defined by Eq. (|96]l . 
The wave function (f>(x) of the bound state is 



( 



D a K 
E 



) 



-l/ct 



4>{x) 



2nh 2 Ea\x\ 




(1 1) (1 - - -) (- -) 



a' an ^2' 2 



(102) 



here the constant C = </>(0) such, that normalization condition 
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dx\(j>(x)\ 2 = 1, (103) 



has to be satisfied. 



5.2 Linear potential field 

Considering a particle in a linear potential field (for example, see jS], page 74), 
the potential function V(x) can be written as: 

ft *>0( f >0) 

OO £ < 0, 

and fractional Schrodinger equation Eq.(f3"T]) becomes 



- D a (7N) a (j){x) + Fx(f)(x) = E(j>(x), Ka<2, x > 0. (105) 

The continuity and bounded conditions of the wave function, let us conclude 
that <p(x) =0, x < 0. Besides, <fi(x) must satisfy the boundary conditions 

f(x) =0, x = 0, 
<p(x) =0, x — > oo. 

Then, wave function <fr n (x) of the quantum state with energy E n , n = 
1,2, 3,. ..is [29] 

tf> n (x) = (107) 



-Hr. 



1 2,2 

a + 1 



pi/ n \ ~^tt f1 3— -J— 1 ( a+2 a 'i 

£jn \ 1 I Ua \ I V a+1 ' a+1 V 2(a + l) ' 2(a + l) - 1 



F'h\(a + l)Fh ' (0 1) ( °+ 2 , , " ) 

(108) 

with the constant A given by 



and the energy spectra E n 



n \ -V(o+i) 



E n = X n Fh^j—^ F ^j , Ka<2, n = 1,2,3,..., (110) 

where A„ are the solutions of the equation [35] 
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H. 



2.2 



-A„ | 



(1 



a+l ' a+1 /' \ 2(a+l) ' 2(cH 
(°> 1 )' (2(a+ 2 l)' 2(a+l)) 



0. 



(Ill) 



When a = 2 Eqs. (|108[ ) and (|1 10[) turn into well-known equations of standard 
quantum mechanics [5], |29j . 

Other solvable physical models of fractional quantum mechanics include ID 
Coulomb potential |29j . a finite square potential well, dynamics in the field of 
ID lattice, penetration through a 8— potential barrier, the Dirac comb [30], the 
bound state problem and penetration through double 6— potential barrier [31J . 



6 Fractional statistical mechanics 
6.1 Density matrix 

In order to develop the fractional statistical mechanics (FSM) let us go in 
Eq. (|40p from imaginary time to "inverse temperature" /3 = l/fc^T, where fcg is 
the Boltzmann's constant and T is the temperature, it — > h/3. Then the par- 
tition function Z is expressed as a trace of the density matrix p L (x, (3\xq) [2], 

M 

Z = dxp L (x,p\x) = 



hp 

j dx J Dx(t) J Dp(r)exp{— ~ J du ^—ip{u)x{u) + H a (p(u),x(u)^ , 

x(0)=x(P)=x 

(H2) 

where the fractional Hamiltonian H a (j>, x) has form ()38|) and p(u), x{u) may 
be considered as paths running along on "imaginary time axis", u = it. The 
exponential expression of Eq. (lll2p is very similar to the fractional canonical 
action (H2"]) . Since it governs the fractional quantum-statistical path integrals it 
may be called the fractional quantum-statistical action or fractional Euclidean 
action, indicated by the superscript (e), 

hp 

S£ ) (P,x)= I du{-ip(u)x(u)+H a (p(u),x(u)}. (113) 



o 

The parameter u is not the true time in any sense. It is just a parameter 
in an expression for the density matrix (see, for instance, [I]). Let us call u 
the "time", leaving the quotation marks to remind us that it is not real time 
(although u does have the dimension of time). Likewise x(u) will be called the 
"coordinate" and p(u) the "momentum". Then Eq. (|112l) may be interpreted 
in following way: Consider all the possible paths by which the system can 
travel between the initial x(0) and final x((3) configurations in the "time" h/3. 
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The fractional density matrix p L is a path integral over all possible paths, the 
contribution from a particular path being the " time" integral of the canonical 
action (I113j) (considered as the functional of the path p(u),x(u) in the phase 
space) divided by 1%. The partition function is derived by integrating over only 
those paths for which initial x(0) and final x{0) configurations are the same and 
after that we integrate over all possible initial (or final) configurations. 



6.1.1 A free particle 

The fractional density matrix p^\x, /3\xq) of a free particle (V = 0) can be 
written as [2j. [14] 



1 

2nh 



(114) 



a\x — xo\ 



tH. 



i.i 



\x-x \ . (1,1/a), (1,1/2) 



h(D a (3)i/« 1 (1,1), (1,1/2) 

where H^j is the Fox H- function (see, [T8 ] -[20 ] ). 

For ID system of space scale the trace of Eq. (lll4p reads 



< fa p w (a .^| a . o) = ^_^_ 7S r ( I ) . 



When a = 2 and D2 = l/2m Eq. (|114p gives the well-known density matrix 
for ID free particle (see Eq.(10-46) of Ref . .1, or Eq.(2-61) of Ref. [32]) 



1/2 



exp 



rn. 2 

^— [X — Xq) 

2h 2 p 



The Fourier representation (p, j3\p') of the fractional density matrix p^' [x, /3\xq 



(115) 



defined by 



Pl\P'P\p') = / dxdx PL\x,P\x )exp{-^(px - p'x )} 



can be rewritten as 



pT(p,f3\p') = 2TThS(p~p')-e-P D ^ a . 

In order to obtain a formula for the fractional partition function in the limit 
of fractional classical mechanics let us study the case when Uj3 is small. It is 
easy to see that the fractional density matrix p L (x, f3\xo) can be written as 



p L (x,f3\x ) = e-P v ^^ I dpexp 



1 

2~h 



'i P -K^~PD a \p\ 
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Then the partition function Z in the limit of classical mechanics becomes 



— oo — oo 

where T(l/a) is the gamma function. 

The partition function Z given by Eq. (|116[ ) is an approximation valid if the 
particles of the system cannot wander very far from their initial positions in the 
"time" h/3. The limit on the distance which the particles can wander before the 
approximation breaks down can be estimated from Eq. (jll4p . We see that if the 
final point differs from the initial point by as mush as 



l/a 



the exponential function of Eq. (|114l) becomes greatly reduced. From this, 
we can infer that intermediate points only on paths which do not contribute 



greatly to the path integral of Eq. (|114[) . Thus, we conclude that if the potential 
V(x) does not alter very much as x moves over this distance, then the fractional 
classical statistical mechanics is valid. 

6.1.2 Motion equation for the density matrix 

The density matrix p^(x, /3\xq) obeys the fractional differential equation [2], |14) 
- 9pL{X g ^ Xo) = -D a (hV x ) a PL (x,P\x Q ) + V(x)p L (x,P\x ) (117) 

or 



- dpL{x ^ Xo) = H a p L (xJ\x Q ), p L (x ) 0\x )=5(x-x ), (118) 

where the fractional Hamiltonian H a is defined by Eq. (TT5l) . The last equation 
can be considered as fractional generalization of the Bloch equation for density 
matrix |33j . 

6.1.3 3D generalization of FSM 

The above developments can be generalized to 3D dimension. It is obviously 
that a free particle density matrix (r,/3|ro) for 3D case has a form 

pf ( r ,£|r ) = -^3 J d 3 P -exp|z P(r 7 r ro) -^ a |p|°|, (119) 
where r, rg and p are 3D vectors. 
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To present the density matrix p L (r, (3\ro) in the terms of the Fox H- function 
we rewrite Eq. (|119l) as 



pf(r,/3|r ) = — -J / d P p S m( Plrb h ^ ) exp {-PD a \p\ a } ■ 

2ir 2 h r — r 7 n 

o 

With help of the identity p<® (r, /3|r ) = -££p£ 0) foj8|0)| B=k _^|,whare 
Pl(x, /3|0) is ID density matrix given by Eq. (|114[) . we find 



|r-r | . (1,1), (1,1/a), (1,1/2) 



h(p a p)V« 1 (1,1), (1,1/2), (2,1) 

(120) 

This is new equation for a free particle fractional density matrix in 3D space. 
The density matrix p L (r, P\tq) obeys the fractional differential equation 



;|r,l) = D a (-h 2 Ar/ 2 p L (r,p\r Q ) + V(r)p L (r,f3\r ), (121) 



dp 



dp L (r,l3\r ) 

d/3 



H a p L (r,l3\v ), p L (r,j3 = 0\r ) = 5(r-r ), (122) 



where 3D fractional Hamiltonian H a is defined by Eq.©. 
Thus, the Eqs. (II 14[) . (| 1 1 (II 22[) are fundamental equations of fractional 
statistical mechanics. 
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